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Abstract 

We present the string dual to finite temperature SU (Nc) = 4 SYM coupled to massless funda- 
mental matter introduced by Nf D7 branes, with Abelian flavor symmetry. The analytic solution 
includes the backreaction of the flavors up to second order in the parameter that weighs the inter- 
nal flavor loops, = {XhNf)/{8Tr'^ N^), \h being the 't Hooft coupling at the temperature of the 
dual Quark-Gluon Plasma. We study the thermodynamics of the system and its departure from 
conformality, which is a second order effect. We then analyze the energy loss of partons moving 
through the plasma, finding that the fundamental degrees of freedom enhance the jet quenching. 
The whole setup is generalized to D3-D7 systems with D3-branes placed at the tip of a generic 
singular Calabi-Yau cone over a five-dimensional Sasaki-Einstein manifold. We finally provide the 
equations for the inclusion of massive flavors in the = 4 SYM plasma. 



1 



1 Introduction 



The "fireballs" experimentally produced at RHIC and the ones that will be produced at the 
LHC contain a non-zero fraction of degrees of freedom in the fundamental representation [1]: 
they are truly Quark-Gluon Plasmas rather than just Gluon Plasmas. While holographic 
methods provide interesting tools to analyze these systems, they have mainly been concerned 
with plasmas without flavors in fundamental representations, or have treated the latter in the 
quenched approximation.^ The simplest and best studied example is the plasma oi H — ^ 
SYM which, unexpectedly, has proven to be not very different from the QCD one. 

This paper is concerned with the study of fiavor effects in the planar, strongly coupled, 
SU{Nc) A/" = 4 SYM plasma beyond the quenched approximation. A solution is presented in 
section 2, dual to A/" = 4 SYM coupled to A^/ 3> 1 massless flavors. The latter are introduced 
by means of D7-branes in the gravity background. The branes are homogeneously smeared 
over the transverse directions of the geometry and so the flavor symmetry group is a product 
of Abehan factors [5, 3].^ The construction is immediately extended to any D3-D7 setup 
with Nc D3-branes at the tip of Calabi-Yau cones over Sasaki-Einstein manifolds. In all the 
cases the zero temperature gauge theories are taken to preserve J\f — 1 supersymmetry in 
4d. 

The solution can be given in closed analytic form up to second order in e^, which is 
essentially XhNf/Nc (precise deflnitions will be given below), where A^^ is the 't Hooft coupling 
at the temperature of the plasma. The validity of the approach requires e/j to be small (see 
section 2.4), which is a common feature of D3-D7 systems. Actually in the particular example 
that we work out in the present paper, realistic inputs yield values of eh ~ 0.24 (see section 
4.1.1 below). Therefore, keeping 0(e|) is already accurate up to the level of a few percent, 
the third and higher order corrections to the solution being tiny. On the other hand, it is 
important to keep the terms since, as we will see, this is the order at which conformality 
breaking (at the quantum level) affects the thermodynamical observables. 

The gauge theories we focus on, in fact, are conformal in the unflavored case and, once we 
couple them with dynamical flavors, have a UV Landau pole. The latter is mapped to the 
blowing up of the dilaton and UV pathologies in the gravity solution at a flnite energy scale 
[7, 8] . In order to make meaningful physical statements, we will take the usual point of view 
of fleld theory: if the IR scale in which we are interested is far below the Landau pole scale, 
it is possible to make well-defined predictions for the IR physics in terms of IR quantities. 
An eventual UV completion affects the results at most in powers of {^ir,uv being the 
IR and UV energy scales). 

The first plasma property we study in section 3 is the thermodynamics. Up to first order 

^Exceptions to this statement can be found in [2, 3, 4]. 

^D3-D7 solutions at zero temperature, where the D7's are locaUzed rather than smeared have been 
discussed in [6]. 
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in XfiNf/Nc our results for the thermodynamic quantities (entropy, energy, free energy, heat 
capacity, speed of sound) confirm the probe computations in [9]. The second order results 
provide new insights into the plasma. In particular, we verify the breaking of conformality 
in the thermodynamic observables at second order, such as the departure of the square of 
the speed of sound from the value 1/3. Since we neglect higher derivative corrections to our 
gravity solutions, the usual relation rj/s — I/Att holds trivially [10], giving automatically the 
value of the shear viscosity. 

An interesting result is provided by the study of the energy loss of partons in the plasma 
in section 4. Using the standard procedures [12, 13, 14, 15, 16] we find that the energy loss is 
enhanced by the fundamental fields. In particular, the holographic value of the jet quenching 
parameter is shifted to higher values with respect to the unflavored plasma, approaching the 
experimental window. 

We finally start considering the inclusion of massive flavors, whose treatment is more 
challenging. We are able to provide in section 5 the equations of motion for the system, but 
we leave their study for future work. A summary of the results is presented in section 6 and 
technical details of calculations are relegated to three appendices. 

2 The D3-D7 plasma solution with massless flavors 

In [7] , the gravity duals to a class of jV = 1 4d quiver gauge theories coupled to massless 
fundamental flavor flclds, were found. The gauge theories describe the low energy dynamics 
at the intersection of A^^ "color" D3-branes and Nf "flavor" D7-braucs. The dual supergravity 
solutions account for the combined backreaction of both stacks of D-branes and thus allow 
an exploration of the non perturbative dynamics of the corresponding gauge theories beyond 
the quenched approximation. 

In the setup the D3-branes are placed at the tip of a Calabi-Yau cone over a 5d Sasaki- 
Einstein manifold X^. In the absence of flavor branes they source a background whose near 
horizon limit is AdS^ x X^: the dual gauge theories are superconformal quivers. Just to pick 
two well known examples: when X5 = the CY manifold is just the 6d Euclidean space 
and the dual fleld theory is A/" = 4 SYM; when X5 = T^'^ the CY manifold is the singular 
conifold and the dual theory is the Klebanov-Witten quiver [17]. Since this will be used 
in the following let us remember that the metric of a 5d Sasaki-Einstein manifold can be 
written as a [/(I) fibration over a 4d Kahler- Einstein manifold: 

dsx, = ds^E + {dT + AkeY , (2.1) 

where r is the fiber coordinate and A^e is the connection one-form whose curvature is related 
to the Kahler form of the KE base: dAxE = ^Jke- For X5 = the KE base is CP^ and 
for X5 = T^'i it is S^xS^. 
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The D7-branes introduce fundamental matter in the dual field theories. They are extended 
along the radial direction of the background and wrap a sub manifold X3 of X5. They are also 
homogeneously smeared over the transverse space [5, 3].^ The smeared distribution is taken 
in such a way that the isometrics of the fibered Kahler-Einstein space arc kept unbroken, 
and allows to write an ansatz where all the unknown functions just depend on a single radial 
coordinate. The D7-brane embedding is taken such that, when the temperature vanishes, 
the whole D3-D7 system preserves M — 1 supersymmetry in 4d.^ As a general feature of 
D3-D7 setups, the dilaton runs and blows up at a certain radial distance, corresponding to 
a UV Landau pole in the dual gauge theory [7] . 

In the following, we arc going to study the finite temperature behavior of these D3-D7 
gauge theories. As a first step, we are going to find the non-extremal generalizations of the 
supergravity solutions found in [7]. 

2.1 Ansatz and equations of motion 

We look for solutions of type IIB supergravity coupled to Nf D7-brane sources. The action 
of the system reads: 

S = Sub + Sfi , (2.2) 
where the active terms of Sub are (we will work in Einstein frame throughout the paper): 



whereas the action for the D7 fiavor branes takes the usual DBI+WZ form: 



(2.3) 



S 



fi 



The gravitational constant and D7-brane tension are, in terms of string parameters: 

1 T7 1 



The metric ansatz we want to focus on is: 



(2.5) 



dsl^ = h--2 [-b dt^ + dxl] + h-^ [b S^F^da^ + S^ds^E + F\dT + Ake)^] ■ (2.6) 



^See [18, 7, 19, 20, 21, 22, 23, 8] for other studies of this kind of construction in the zero temperature 
case. 

^When X5 = S^, a. single stack of parallel D7-branes preserves = 2. However, the smeared distributions 
considered in [7] and in the present paper just preserve J\f = 1. See appendix C for more details on the 
D7-brane embeddings for this case. 
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The functions h, b, S, F (as well as the dilaton $) depend on the radial variable a; b — 1 
in the extremal (i.e. zero-temperature) case. In our conventions, S, F have dimensions of 
length, b, h are dimensionless and a has dimension length"^. 
The non-trivial RR field strengths are: 

i^(5) = Qc (1 + *)e{X^) , = Qf {dr + Ake) , {dF^^l) = 2 Qf Jke) , (2.7) 

where £{Xr^) is the volume clement of the internal space and Q^, Qf are proportional to the 
number of colors and flavors: 

{2nr 9s a- ' ' Vol{X,)g. " ^^'^^ 

The first relation is just the usual quantization condition of the D3-brane charge, while the 
derivation of the second is detailed in appendix A.^ Vol{X-i) is the volume of the submanifold 
(X3 C X5) wrapped by any D7-brane.^ In the X5 = case, Vol^X^) = 27r^. 

The fact that the flavors are massless is encoded in the independence of on a. See 
[7, 20] for discussions on this issue and section 5 for the case of massive flavors when X^ — S^. 

All the functions we need to compute depend on a single coordinate a, and it is possible 
to describe the system in terms of a one-dimensional effective action. By directly inserting 
the ansatz in the action (2.2), we find: 

^ 24, J + ^^^^ + 2^^^^*^ + 

(d^b) f (d^h) (d^S) {d.F)\ 

+ 24bF^S^-4bF^S^-^Ql-^-^Qy%S'' -AQfe^'bS^'F'^^ . (2.9) 

In (2.9) Vi^3 denotes the (infinite) integral over the Minkowski coordinates. The second 
derivatives coming from the Ricci scalar have been integrated by parts and, as is customary, 
only the angular part of F(5) is inserted in the F^^^ term (otherwise the Qc would not enter the 
effective action since, on-shell, F^^-j = due to the self-duality condition). Notice also that 
the WZ term docs not enter (2.9) because it does not depend on the metric or the dilaton; 
its effect has been taken into account via the expression for (we refer the reader to [7] 



'^Wc have defined Nj- as the number of flavor branes. The precise way in which they introduce fundamental 
degrees of freedom and how they couple to adjoints and bifundamentals depends on the particular theory. For 
instance, in the = case, there are 2Nf supermultiplets (in A/" = 1 language) but the flavor symmetry 
is just U{l)^f , due to the coupling between adjoint and fundamental superfields in the superpotential. 

^Notice that, since the different D7-branes take different positions in the internal space due to the smear- 
ing, the particular submanifold X3 spanned by each brane is difi^erent. However, since the family of relevant 
X^'s are related by the internal isometrics, they all have the same volume and thus Vol{X3) is a well defined 
quantity. 
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for extensive explanations). The term proportional to Qf comes from the DBI contribution 
in (2.4); more detailed considerations regarding this last term are relegated to appendix A. 
The equations of motion stemming from the effective action (2.9) are: 

dlilogb) = 0, 

dli\ogh) = -Ql^, 

dlilogS) = -2hF^S^ + QhF^S'' - Qf e%F^ 5^ , 

dl{\ogF) = AhF^S^-'^e''HS\ 

dl^ = Q}e''^bS^ + AQfbe^S''F\ (2.10) 

It is straightforward to check that these equations solve the full set of Einstein equations 
provided the following "zero-energy" constraint is also satisfied: 

= 4(^,i,(^,s(^^-l(...)^. 

(d^b) f (d.h) (d^S) {d„F)\ 

- 24bF^S'^ + 4bF^S^ + -Ql^ + -Qje^%S^ +4Qfe''bS'^F\ (2.11) 

This constraint can be thought of as the aa component of the Einstein equations or, alter- 
natively, as the Gauss law from the gauge fixing of g^rcr in the ansatz (2.6). By differentiating 
(2.11) and using (2.10), one finds an identity, which shows that the system is not ovcrdctcr- 
mined. The equations are valid for any Kahler-Einstein base space since only the properties 
R^ff — QQab^ and dAxE — '^Jke are needed when checking the Einstein equations. 



2.2 The supersymmetric solution: second order expansion 

The supersymmetric (zero temperature) solution of (2.10), (2.11) was found in [7]. It corre- 
sponds to having 6 = 1 and it is dual to A/" = 1 quiver gauge theories with flavors. In the 
present notation, the BPS first order equations read: 

d,^ ^QfS^e^, d,F ^S'F(^3-2^-^e''^ . (2.12) 

It can readily be shown that these equations, together with 6=1, arc sufficient conditions 
for (2.10), (2.11) to hold. In a different, dimensionless, radial coordinate dp = S^da, the 
equations for S, F, $ can be explicitly integrated, giving [7] : 

S = a'^e" (l + e*(J + p*-p) 
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$ = $*-log(l + e*(p*-p)) , 



^ = _g^«'-2e-4p/i + ,^(_ + ^^_^)j ^ (2.13) 

where for later convenience a scale has been introduced and is the value of the dilaton 
at that scale. We have also inserted powers of a' (which enters as an integration constant 
of (2.12)) in order to give appropriate dimensions. Notice that the solution is defined for 
P < Plp where pip = p* + e;^^ is the point at which the dilaton blows up. We have 
chosen to keep the differential equation for h, even if it can be solved in terms of incomplete 
gamma- functions. The parameter e* = Qf^^* has been introduced. It has to be small 
for the solution to be valid in a large energy range (see section 2.4) and it will be used as 
an expansion parameter. Defining A* as the 't Hooft coupling^ at the p* scale, e* can be 
expressed in terms of physical quantities by using (2.8) as: 

- -A.t , (2.14) 



16n Vol{X5) 

and, in particular, 6^^X5=8^) = s^^*^- 

We have set to the integration constant Ci of [7] for the sake of IR regularity. Even if 

for Ci = 0, this backreacted geometry still presents an IR singularity (much milder than in 
the Ci 7^ cases), it is useful to think of the ci = solution as the massless limit of a family 
of IR regular solutions where the IR singularity is removed by non-zero quark masses [20]. 
When we go to finite temperature, the singularity will be hidden behind an event horizon. 

In comparing to the unflavored and to the finite temperature solutions, it will be useful 
to employ an r coordinate which we define by requiring that h takes the simple and familiar 
form: 

h-^ = 1^ _ 1^ (2-)V^ .215) 

We can expand ^ from (2.13) and integrate order by order in e*. Since (2.15) gives h 
explicitly in terms of r, this yields an expression for r(p). Let us fix the additive integration 
constant in h such that r(p^) = = \fa'eP* ? Then: 

^1 + 1^ (^,4p-4p. _ 1 + i2(p, _ p)) + -M- (e«P-8p. + ^e'^-^^^ (3 + 4(p. - p)) 



For the (flavored) A/" = 4 SU{Nc) theory, the gauge coupUng is — ^ti" ffsC* (note the choice for the 
numerical prefactor, which sometimes is taken to be 27r), and thus A* = 47r ^s-^cC** ■ For quiver theories 
that correspond to difii'erent X5 geometries, the gauge groups are of the form SU{Nc)'^. Let us generaUze a 
relation from the orbifold constructions AirgyM ^ = (g^e*)"^ [17, 24], and consider all the gauge couplings 
9YMA to be equal. Then in ggNce"^, strictly speaking, gives the 't Hooft coupHng at each node of the quiver, 
divided by n. With an abuse of language we will simply refer to it as the 't Hooft coupling. 

^This is different from the choice of integration constant adopted in [8], which required h to vanish at the 
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-(19 - 24(p. - p) + 144(p, - + 0{el)\ . (2.16) 
We can obtain now F{r), S{r), $(r) as expansions up to second order: 

24^ 3rr 1152 V 9 9r8 ^rj J ^ *. 

2 / 4 

r ez I ^ r r nr 



-So = r 

$0 = $. + e.log- + ^ (l-^ + 121og- + 361og^-) +0(e^) , (2.17) 
where the subscript means that the solutions are dual to the D3-D7 theories at T = 0. 

2.3 The non-extremal solutions 

We now look for non-extremal solutions of (2.10), (2.11), which would provide for a dual 
description to the finite temperature regime of our D3-D7 gauge theories. Such solutions 
are required to be regular at the horizon and to tend to the supersymmetric (T = 0) ones 
at energy scales much higher than the temperature. Concretely, we will require that the 
geometries coincide with the T = solutions in the extremal limit and that F, S*, $ coincide 
with those in (2.17) when evaluated at r = r*. This uniquely fixes the order by order 
expansion of the non-supersymmetric solution. 

The equations for h and 6 in (2.10) are decoupled from the rest and are solved, in terms 
of an integration constant r^, by: 



h^e<'' , h^^(l-e<''). (2.18) 

4r? 



h 



where a G (— oo,0). We define the r coordinate such that the expression for h (2.15) still 
holds: 



4 

The extremal limit corresponds to sending the horizon radius r/j to zero. The metric reads: 



e^^t'^ = 1 - _ (2.19) 



dslo--^A^-'4)dt' + ^2d4+ ^'K^' +R'S'dsj,^ + R'F^{dT+A^^r , (2.20) 

K r K r (^ — ^) 

where: 

5 = - , F= - . (2.21) 

r r 



point where the dilaton diverges. However, the difference between both integration constants is a quantity 
±_ 

suppressed as e and its influence on the IR physics is of this order and therefore neghgible with respect 
to the terms kept in the Taylor expansion. 
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We still have to solve for F,S,^. A straightforward computation from (2.10), (2.11) leads 
to the differential equations for F, S, $, in terms of the coordinate r. It is easy to check 
that the AdS^ black hole solution F = S = 1, ^ = const is recovered in the flavorless limit 
= 0. Its deformation, expanded up to second order in reads:^ 

F = 1 ^*ni ^^'~^M l (l- 94 2r^-r^ 5 (2r^-r^)^ 

24^ ^ 6rt -3rt'^ 1152 V 9 2rt - 9 {2rt - r^)^^ 



9 {2rt - 



-481og(-) +0(e 



_ 2r^ - rj 106 2r^ 5 (2r^ - rif 

24^ 6rt-3rl' 1152 \ 9 2rt - 9 {2rt - rl)"^ 



8/4 4\ 



where Li2{u) = Yl'^=i ^ is a polylogarithmic function. The expressions in (2.20), (2.22) 
are the central results of this paper. Together with (2.7) they provide the full perturbative 
solution. We are now going to discuss their regime of validity and, in the following sections, 
their implications for the physics of the dual flavored plasmas. 



2.4 Hierarchy of scales and regime of validity 

In order for the set-up to be physically meaningful, there must exist a hierarchy of scales. 
In terms of the r radial coordinate: 



< r* < r„ < rLP ■ (2.23) 

The quantity Vh sets the scale of the plasma temperature ^ ~ Ajr ~ T, which is the scale at 
which we want to analyze the physics. The point r^p is where the dilaton diverges, signaling 
a Landau pole in the dual theory. However, as discussed in [8], the string solution starts 
presenting pathologies at a lower scale, which we denote by r^. At this scale, which is fairly 
close to r^p, the holographic a- function is singular and the utility of the solution for r > Va 

^This procedure is analogous to the one adopted in [25]. 

^°Bcing strictly precise, there should be an extra scale < <C r/j proportional to the quark masses 
if one wants to avoid the IR singularity of the zero temperature solution. The solutions presented in this 
section can be seen as the leading ones in ^ ~ If one wants to use the results for phenomenological 
estimates of the QCD plasma, this should be a good approximation since m„, rrid -C T. Section 5 deals with 
the backreaction of massive quarks. 
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is doubtful. Finally, r* sets an (arbitrary) UV cutoff scale ^ ~ ^uv- The solution (2.22) 
will only be used for r < r*. In a Wilsonian sense of a renormalization group flow, the UV 
details should not affect the IR physical predictions. This feature is reflected in that physical 
quantities do not depend (up to suppressed contributions) on r^, or functions evaluated at 
that point, but only on IR parameters. Even if the precise value of r* is arbitrary, we have to 
make sure that it is possible to choose it such that it is well above the IR scale (so that the 
UV completion only has neghgible effects on the IR physics) and well below the pathological 
fa-ifLP scales (so that the solution we use is meaningful and the expansions do not break 
down). To this we turn now. 

Let us start by computing the hierarchy between r* and tlp- Since at we can approxi- 
mate the solution by the supersymmetric one, we can read the position of the Landau pole 
from (2.13). If we insert the approximate relation between radial coordinates r a; \fa'e'' ^ we 
find: 

— e~^^ < 1 , (2.24) 

as long as e* ^ 1. 

Moreover, one has to make sure that the Taylor expansions 
Tfe < r < r*. This of course requires e* -C 1, but also that e* 



2.22) are valid in the region 



log 



<C 1 (notice that the 



absolute value of the logarithm can be big because r/, ^ r*). This means that ^ e . On 
the other hand, when we compute physical quantities in the upcoming sections, we always 



neglect quantities suppressed as powers of — ~ t^-^^ This is the order of magnitude of the 



Al 



corrections due to the eventual UV completion of the theory at r\. One has to make sure 
that the corrections in e* we are keeping are much larger than the neglected ones, namely 
e* ^ — . In summary, we have the following hierarchy of parameters (in the following, in 
order to avoid overly messy expressions, we insert the value of e* for the = case, 
remembering that for a generic X5, its value is given by (2.14)): 



~ e ^ « ^ - « e. ~ « 1 . (2.25) 



As long as ~ <^ 1, there always exists a range of r* such that this inequality is 

satisfied. Since we will focus on the IR physics of the plasmas, at the scale set by their 
temperature, the actual physical constraint on the parameters will be ^ 1 , which we 

have written in terms of the coupling at the scale of the horizon, = A*(l + 0(e*)), see 
section 3 below. 

On top of this, we have to make sure that the SUGRA+DBI+WZ action we are using 
is valid. As usual, the suppression of closed string loops requires Nc ^ 1 whereas the 



^^For the thermodynamical quantities, we always neglect terms suppressed, at least, as (j^^ ; while for 
the jet quenching parameter the neglected terms are of order — . 
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suppression of ct'-corrections is guaranteed by A/j ^ 1. We have written the D7 worldvolume 
contribution to the action as a sum of Nf single brane contributions. This is justified if the 
typical energy of a string connecting two different branes is large (in a' units). Since the 

branes are distributed on a space whose size is controlled by i? ~ A/^v^, we again need 
A/i ^ 1. The smearing approximation will be good if the distribution of D7-branes on the 
transverse space is dense, i.e. Nf ^ 1. The discussion up to now is summarized in the 
following validity regime: 

TV, > 1 , A^ > 1 , TV/ > 1 , en^ -\-f- < 1 . (2.26) 

Finally, we want to find the regime of parameters in which the fiavor corrections are not only 
valid but arc also the leading ones. With this aim, we ought to demand that a'-corrcctions 

_ 3 

to the supcrgravity action (which typically scale as A^^ ^ due to first string corrections of the 
type a'^R^) are smaller than the fiavor ones, controlled by eh, namely:^^ 

A^^ « e, . (2.27) 

Demanding that corrections to the D7-branes contributions (for instance curvature correc- 
tions to the worldvolume action itself or corrections produced by possible modifications of 
the brane embeddings due to curvature corrections to the background metric) are subleading 
does not impose any further restriction. The reason is that their contribution is typically of 
order e^^A^"^ for some c > which is always subleading with respect to €h as long as (2.26) is 
satisfied. 

Let Tis conclude this section with some comments on the stability of our perturbative 
non-cxtrcmal solutions. A possible way to check for the latter is to consider worldvolume 
fiuctuations of a D7-brane in the setup. If, as in our cases, the brane corresponds to massless 
flavors, the related quasi-normal modes on the unflavored background all have frequencies 
with a negative imaginary part of the order of the temperature, signaling stability [26], [9]. 
This result cannot be changed in the flavored case when a perturbative expansion in eh is 
done. Thus, in our regime of approximations, stability with respect to those fluctuations is 
guaranteed. 



3 Thermodynamics of the solution 

In the previous section we have found a family of perturbative non extremal solutions which 
are regular at the horizon. They are dual to finite temperature fiavored deformations of 
conformal theories, including A/" = 4 SYM. The solution (2.22) is written in terms of the 

""^^Notice that if one wants the second order flavor corrections that will be computed to also dominate over 

_ 3 

the curvature ones, the more restrictive condition ^ ^® needed. 
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parameters e*,r*, which are defined at the UV scale. As already stressed, the physical 
quantities must be expressed in terms of IR parameters. We thus define: 

_ KVoljX,) Nf 

where the subscript h means that the quantities are evaluated at the horizon. Thus, Xh is 
naturally identified with the coupling at the scale of the plasma temperature. We therefore 
have: 



= = e* + e* log h 0(ej . 



(3.2) 



It is important to notice that, if we tune the temperature while keeping the UV parameters 
fixed (namely, we change the temperature without changing the theory), then depends 
on the temperature. Since, as we will see below, Vh is proportional to the temperature (at 
leading order), we have: 

(3.3) 



and T{dXh/dT) = ehXh at leading order. These relations refiect the running of the gauge 
coupling induced by the dynamical fiavors. 

We now ask how the BH temperature and entropy density are related to the parameters 
of the solution. Let us start by computing the temperature requiring regularity of the 
euclideanized metric, by identifying the temperature with the inverse of the period of the 
euclideanized time. A simple computation yields: 



2r. 



(3.4) 



where we have inserted the values of F,S at the horizon, which can be read from (2.22) 
(neglecting terms suppressed as powers of -|): 



, eh 17 

1 - — H 

24 1152 



el + 0{el) , ~Sh 



The entropy density s is proportional to Ag, the volume at the horizon of the eight dimen- 
sional part of the space orthogonal to the r plane (where t is the Euclidean time), divided 
by the infinite constant volume of the 3d space directions V^. Prom the general form of the 
metric (2.20) we get that: 



27r ^8 rlR^SlFhVol{X^) 



which in terms of the temperature reads: 



2Vol{X^) "n^R^^ 



1 19 



2Vol{X5 



(3.6) 



(3.7) 
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As for the other thermodynamic quantities which will follow, the leading term of this formula 

is the well-known unfiavored result. The 0{eh) term was already calculated in [9] with the 
probe brane technique, in the X5 = case. We re-obtain the result in a quite standard 
way, by computing the increase of the horizon area produced by the flavor branes. This can 
be considered as a crosscheck of the validity of the whole construction. Finally, the order 
term is a previously unknown contribution. 

The expression we have found in (3.7), shows that the relation s — 27r'^aFT{T)T^, where 
apT is the holographic a-charge of the T — theory, holds also in the flavored case, al- 
though only up to first order in eh- Let us show this exphcitly. The general formula 
for the holographic a-function of a flavored T = SQCD-like theory dual to the D3-D7 

background is given hy a pT = 

27{8Tr^gy^)-^l32H2{drH)-^ [27] where we have defined 
H = R^r^VoliX^fS^F^ and /3 = KS^Fq [28]- Inserting (2.17), we get: 



apT — 



^Vol{X^) 



1 + + e 
2 



1 1 r 5 r 
6^2 972^ 



+ 0{el) 



(3.^ 



This formula shows that the holographic central charge is constant up to first order in e*. 
However, it starts running at second order, as a signal of the broken conformal invariance. 
Defining apT{T) — apxirh), we see that the claimed relation with the entropy holds only at 
first order. 

The ADM energy of the solution can be straightforwardly computed (see appendix B). It 
yields the energy density of the plasma and, thus, it allows us to study the full thermody- 
namics: 



E 



e — 



ADM 



-NlT^ 



1 + leniT) + \eu{Tf + 0(e,(r)=^) 



(3.9) 



Again, terms suppressed as powers of ^ have been neglected. Moreover, since in the following 
derivatives with respect to T are going to be taken, we find it convenient to make explicit 
that e/j depends on T (see (3.3)). From the relation above we get immediately the heat 
capacity (density): 



2Vol{X5) " 



l + ^e.(T) + lie.(T)2 + 0(Q.(m 



(3.10) 



The free energy density, and so (minus) the pressure, reads: 



8 Vol{X5 



1 + \eH{T) + \en{Tf + 0{en{Tf) 



(3.11) 



Notice that, consistently, this satisfies the relation s = OtP (where it is crucial to take (3.3) 
into account). This result is confirmed by the direct computation of F, which is relegated to 
appendix B. Comparing (3.11) to (3.8), we find that p = \apT{T)T'^ up to second order. 
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The speed of sound Vg is: 

e 1 r 1 _ „ 1 

(3.12) 



9 S 1 

Cv o 



le,{Tr + 0{e,{Tf) 



Note that the correction to tlic speed of sound only appears at second order. Instead, when 
quarks are massive, they break conformal symmetry at tree level and the correction to the 
speed of sound is of first order in [9]. It is also interesting to notice the sign of the 
correction, which is consistent with the upper bound Vg<^ conjectured in [29] for a certain 
class of theories (see also [30]). In [29], a heuristic motivation of the bound was given on 
the basis of asymptotic freedom, which it is not the case for the present theory. As far as 
wc know, (3.12) is the first explicit check of this bound in the gravity dual of a theory with 
positive /3-function. 

Analogous to the running of gft (eq (3.8)), the deviation from conformality in (3.12) is a 
second order effect. In fact, the solution provides a measure of the breaking of conformality 
at second order from the interaction measure: 

MTf. (3.13) 



16^0/(^5) 

The determination of this as well as other thermodynamic quantities on the lattice for QCD 
can be found in [31]. Of course, the running of observables in QCD and in our D3-D7 models 
are opposite. 

Since we have not introduced into the action terms with higher derivatives of the metric 

(such as curvature squared terms), the usual theorems apply and the fundamental matter 
does not affect the result 7 = 4^ [10]. Wc can thus read the value of the shear viscosity 
trivially by dividing (3.7) by Arr. Again, while the first order value was already calculated in 
[33], the second order result is new. 

Concerning the bulk viscosity, if our thermal gauge theories saturate order by order the 
bound proposed in [34]: 

|>2(l-«j). (3.14) 

then C would still be zero up to first order. On the other hand, at second order we would 
obtain a non-trivial result :^'^ 

C = j2Vd{X,) ^'^' + 0(6,(T)3)] . (3.15) 

Adapting the reasoning in [35] to our case, one can find indications that the bound is indeed 
saturated, provided some (reasonable) assumptions are fulfilled. It would be very interesting 



^'^In some sense, our solutions could represent a realization of the "magnetic component" of the QGP as 

described in [32]. 
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The fact that the bulk viscosity should be seen as a second order effect was already observed in 
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to perform the precise holographic calculation of C,/r] to check whether this is actually the 
case. 

4 Energy loss of partons in D3-D7 plasmas 

Let us now consider the energy loss of partons in the quark-gluon plasmas dual to the 
non-extremal D3-D7 solutions we have found in section 2.3. 

The holographic study of energy loss has turned out to be quite relevant from a phe- 
nomcnological point of view. This is because the real-world QCD quark-gluon plasma, 
whose properties are being studied at RHIC and will be investigated at the LHC, appears 
to be strongly coupled [36]. Moreover, the phenomenon of jet quenching observed at RHIC 
demands for a very efficient mechanism of energy loss. There are two main ways to account 
for this phenomenon in the stringy picture. 

On one hand, it is natural to try to model jet quenching as a result of the bremsstrahlung 
that occurs through the strong interactions of the parton probe with the quarks and gluons 
in the plasma. In perturbation theory, this mechanism is effectively captured by a transport 
coefficient termed g, the jet-quenching parameter [11]. At very high energy, and using 
the eikonal approximation, the authors of [12] found a non-perturbative prescription for 
calculating q as the coefficient of in an almost hght-like Wilson loop with dimensions 
L~ ^ L. Being non-perturbative it can be implemented in a string theoretic framework and 
this was done for the first time in [12]. Following this prescription,^® we will compute q for 
our backreacted background in section 4.1. 

On the other hand, at strong coupling the energy loss of a probe passing through a medium 
can be modeled entirely within a string theoretic framework [13, 14]. A parton of velocity v 
is described by a macroscopic string attached to a probe flavor brane. The string is dragged 
by a constant force / which keeps the velocity fixed. The drag force transfers energy and 
momentum to the parton, which are lost in the plasma at a constant rate. The drag coefficient 

which measures the energy loss, is then calculated from the equation / = //p, where p 
is the parton momentum. In section 4.2, we will deal with this drag force in the flavored 
background. 

It is worth noticing that the computations we will make with the Nambu-Goto action 
receive stringy corrections of order A~2 [39]. If we want them to be subleading with respect 
to the flavor corrections, we need to require A~2 <^ e/i, which is more restrictive than (2.27). 

^^See also [37]. 

^^In [38] a different approach to the same problem has been recently proposed. 
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4.1 The jet quenching parameter 



Let us compute the jet quenching parameter q following the prescription of [12]. Taking the 
generic formula in [16]^^ (and cutting the integral at r*), we can write: 



-1 / / -* 
— na 6 2 



9 XX y/ 9 XX + gtt 



ar — 7i — e 2 



'I 

J n 



e 2 — dr . 

rh 



(4.1) 



The dilaton enters the formula because we are considering the Einstein frame metric. We 



have extracted e" 2 because we want to factor out the physical IR parameter — 
e^V^TT 9sNc — e-2^:3!. V^Efp^ can now insert (2.22) and write everything in terms 
of eh rather than e*. Finally, performing a change of variable ^ = we can write the 
inverse of the jet quenching parameter as: 



hVA/j 7r2 



(4.2) 



The integrals that appear at each order are: 

1 - 



h 



i: 

r°° 1 

Ji 8 



- 4 log ^ 



dg = 
dg 



r(!) ' 

V^r(f) l-TT 

r(!) 



87^"^" r(|) 8 

19 - 8 log ^ - 48 log^ ^ - 12 Li2(l - g'^] 



384y? 



dg ~ 



v^r(l) 1 

r(f) 384 



19-48C-7r(2 + 37r)-8 4F3 ( 1, 1, 1, -; -, 2, 2; 1 



3 7 



(4.3) 



where C ~ 0.91597 is the Catalan constant. The upper limit of the integral has been taken to 
infinity instead of = ^ since all the integrands are of order g~^ at large g and, therefore, 
J°° — J^* — 0(g~^) — O(^), which, as usual, we disregard. The jet quenching parameter in 
terms of gauge theory quantities reads: 



rp3 



0M^r(f) 

where we have introduced a constant 7: 



l + -{2 + 7r)en + ^el + 0{el) 



(4.4) 



11 TT 37r2 1 1 ^ 

^= 96 + 48 + 128 + 8^+48 ^"^^ 



1,1,1,|^,2,2;1 



0.5565 . 



(4.5) 



^''We took into account a different factor of between the definition of [12] and of [16]. 
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4.1.1 Possible implications for phenomenology 



We can now discuss some physics coming from (4.4). In [40], it was shown that for a class 
of theories with only adjoint (and bifundamental) fields with gravity duals, one can write: 



r(!) 



(4.6) 



It was argued that deviations from this formula could come from having fundamental fields 
or from non-conformality of the theory. In our case we actually deal with both of these 
effects. As far as we know, this is the first set-up in which one can directly test the effect of 
fundamental fields on the jet quenching parameter in a framework completely under control. 
The expression (4.6) is modified to: 



TT 

8 



TT 



eh 



(4.7) 



Thus, from the point of view of this formula, q increases if Nf > 0. At first sight, one could 
think that this enhancement of the jet quenching parameter is a trivial effect due to the 
increase of degrees of freedom of the plasma because of the addition of new fields. This 
is a naive conclusion since, as remarked in [12], the expression (4.6) shows that, contrary 
to previous expectation, q is not linked to the number of degrees of freedom. However, it 
is instructive to make the following gcdankcn experiment :^^ imagine we have a flavor- /ess 
plasma at temperature T, entropy density s and with A^^^i colors. One can think of replacing 
some of the adjoints by fundamentals, while keeping the total number of degrees of freedom 
(namely, the entropy density s) fixed. We will also maintain fixed the temperature T and the 
coupling constant = gyM/i^^) scale set by T}'^ We thus compare the theory with 

A^c,i colors and no flavors to a theory with Nc^2 (< A^^i) colors and Nf flavors. Concretely, 
keeping T, s fixed requires Nc^2 = N^^i)- — \eh — ^^D? where we have used (3.7). For the 
flavored theory, but written in terms of Nc^i (the corresponding number of colors of the 
unfiavorcd theory), we would have: 



1 9ym 



, l + TT , 25 

1 + —eh + (7 - 



TT 



8 



384 64 



i4 + o(6D 



(4.8) 



The positive sign of the correction terms shows that the fiavors enhance the jet quenching. 
This conclusion is confirmed in the alternative comparison scheme proposed in [41], in which 
the energy density and the force between external quarks are kept fixed, while T and Xh are 



^^We thank David Matoos for suggesting this comparison. 

^^Fixing A/i instead of as produces the same result. Fixing Nc and et gives instead the opposite behavior, 
but it seems a pathological choice, since it does not include the unflavored limit en 0. 
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varied. ^° The fact that the jet quenching is enhanced in all of these comparison schemes (we 
found no reasonable counter-examples) seems to indicate that the result is quite robust. 

Interestingly, the enhancement of the jet quenching by flavors was already observed in [4] 
in the context of a non-critical five- dimensional string dual of a flavored plasma. Of course, 
in that model the string corrections are not under control, while in the present case the result 
is completely trustworthy. 

Let us close this section by making a numerical estimate of the jet quenching parameter, 
by inserting quantities relevant at RHIC, even if it is not clear which is the best way to 
extrapolate the results from the theories here discussed when giving estimates for QCD. As 
a first step, in order to compare with the unflavored jV = 4 SYM result in [12], let us pick 
VollX^) = Vol{S^) — TT^ and extrapolate our result to the realistic regime (which is not 
included in the regime of validity of any holographic model) where ~ 1/2 and A^^c = 3, 
i.e. A/j ~ Gtt. Then wc would have eh ~ ~ 0-24 for Nf = 3. From (4.4) wc read that 

the correction with respect to the J\f = 4 SYM result corresponds to an increase of about 
20% for Nf = 3. As an example, at T = 300 MeV we would get q ~ 5.3 (Gev)Vfm, to 
be compared with the value q ~ 4.5 (Gev)^/fm of the unflavored plasma [12]. The flavored 
result would be right in the ballpark of RHIC values, i.e. g ~ 5 — 15 (Gev)^/fm. 

This comparison with QCD involves two theories with a different number of degrees of 
freedom. Should we compare the two theories at flxed temperature, entropy density and as, 
as discussed above, we would get a smaller q. Let us take sqcd/T^ ~ 17.5 at T = 300 McV 
from [31] (figure 10). Then, using (3.7), and inserting Nf = 3 and eh = wc get a smaller 
coupling Xh ~ 11.1. This would imply, from (4.4), that q ~ 4.1 (Gcv)^/fm. In the same 
scheme, in the flavorless Af — A SYM limit we would get q ~ 3.5 (Gev)^/fm. 

Let us note that the interaction measure in the latter scheme is rather small, {e — 3p) /T^ — 
A^^7r^e^/16 ~ 0.11. In lattice QCD it is also known to be small at temperatures above 2Tc ~ 
400 McV, where the theory is nearly conformal, but not that much. For example, from figure 
4 in [31], it is around 1 for temperatures around 500 MeV. Of course, in QCD {e — 3p)/T'^ 
has a large contribution also at the flavorless level, being proportional to the trace anomaly. 

4.2 The drag force 

Let us now consider a heavy quark moving through our D3-D7 plasmas and compute the drag 
force it experiences, following the general procedure described in [13, 14, 15]. We consider 
a sim]:)le string configuration representing a test quark moving in a given spatial direction 

■^"To be more specific and working at first order in e^: Nc is kept fixed, so from the matcliing of e one gets 
T2 = Ti{l — eh/8), so 32/91 = \^ ^h,2 / ^h,i[^ -|- (tt — l)€h/8]. The couphng A^, is adjusted in such a way that 
the force between two external quarks at the screening length of the first plasma, aqq = 3L^ iF'(Lc,i)/4, is 
kept fixed [41]. By numerically computing the potential V{L) from the Wilson loop in the two plasmas and 
equating the two forces aqq, one straightforwardly finds that Aft,, 2 > ^h,i and so q2 > qi- 
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x:t — T, r — a, X — x{a,T). In particular, we will just discuss the stationary string 
configuration (an open string with an extremum attached to a probe D7-brane at r = r*) 
corresponding to a quark which moves at constant velocity, such that the energy loss due to 
friction with the medium is compensated by an external force. This is achieved by setting 
x{a, f) = r{a) + V f. We have to analyze the Nambu-Goto action in the background metric 
dslfj. — e^/'^ds^ with ds^ given in (2.20). We can use general results from [15] and define 
C as the constant determined from the equation gxx{fc)gtt{fc) + = with the point Tc 
given by gtt{rc) + gxxifcjv"^ — 0, namely Tc — r/i(l — v^)'"^. Then, the rate of momentum 
transferred to the medium is given by [15]: 

-C^-- — 2 = -n Mkin —F , (4.9) 



dt 27ra' 27: a' VT 



V 



where we have introduced notation from [13]: a friction parameter /x such that = — /ip 
and a kinematical mass Mkin such that p = Mkin^j^^- From (4.9), using (2.15), (2.22), 
(3.2), (3.4), we find: 



kin 



l + l(2-log(l-^;2))e^+ (4.10) 



2 ^Vol{X,) 

+ ^ [44 - 201og(l - v') + 91og2(l - v') + 12L^2(t;')] el + O(e^) 

As in section 4.1, the energy loss (at fixed v) is enhanced by the presence of fundamental 
matter, also in the different comparison schemes described in section 4.1.1. The quantity 
fi Mkin grows when increasing the velocity. From (4.10), formally, it would diverge as v —>■ 1. 
However, (4.10) is not applicable in that limit since we have to require ehlog(l — v^) -C 1 for 
the expansions to be valid. 



5 Deforming AdS^ x with massive flavors 

In this section, we will write down the generalization of equations (2.10), (2.11) for the case 
in which the backreacting flavors are massive (all with the same modulus of the mass). The 
search for solutions and the study of the physics is left for future work. Some technical 
details associated to this section are relegated to appendix C. 

The extra complication we face when the quarks are massive is that the D7-brane embed- 
dings are in this case non-trivial and, in fact, have to be studied numerically (when T > 0, 
even if backreaction is not taken into account). Moreover, this case cannot be studied with 
the same generality as the massless one, since the mentioned embedding equations do depend 
on the metric of the Kahler-Einstein space. We will focus in the case where the set-up is a 
flavor deformation of the AdS^ x black hole. Thus, the Kahler-Einstein space is CP^ and 
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the ten dimensional metric we consider is (2.6) with: 

dslp2 = 7^X^ + 7 cos^ ^{de^ + sin^ ed(p^) + \ cos^ | sin^ Md^) + cos Qd^ff , 
1 y 

Acp2 = -cos^ -((i-^ + cos6'd(/7) . (5.1) 

The range of the angles is < x, ^ < tt, < </?, r < 27r, < ^ < 47r. The F(5) RR field 
strength takes the same form as in the massless case, but picks up a dependence on the 
radial coordinate: 

i^(5) = Qc (1 + *)eiS') , = QfPia) (dr + Acp^) , (5.2) 

where p{a) is a function which depends on the branc cmbcddings. It has to vanish at scales 
smaller than the quark masses and to asymptote to 1 at energy scales much larger than the 
quark masses. 

First of all, we have to write down the equation that determines the embedding of a D7- 
brane in this background. Let us take as a "fiducial brane" one wrapping 9,ip,il!, with x a 
function of a and situated at fixed r. The rest of the family of embeddings needed for the 
smearing are obtained from the fiducial one by acting with the symmetries of the internal 
space, see appendix C for a discussion. The DBI action for this fiducial brane reads: 



Snsi ^--jdxe bsmeS F cos [—) ^ cos- [—) + -sm (— . , ^^^^ , 

(5.3) 

where Xwv is the function of cr that determines the brane embedding. The 8d integral is taken 
along the Minkowski directions, 6, (/?, ip and a. There is also a WZ term due to the coupling 
to the background C(8). The C(8) potential can be computed from (iC(8) = e^*(*F(i)). The 
WZ piece of the action is: 

Swz = ^^Q/ J d'xsmep(a)be'^S'cos' (^) . (5.4) 
The equation of motion for Xwv follows from the action Sy^y — Sdbi + Swz' 
= la,(^e*cos^^|i(a,x«,.))+ (5.5) 
+ e*65«F2cos2^sin^(^3H,S2 + cos2^(l-|J)| + g,e*|^cos^p(a^^ . 
In order to abbreviate the notation, we have introduced the quantities: 



/ 9 / Xwv\ , S'- . 2 fXwv\ ^ /1 I {9aX'WVJ n\ 

= a/cos^ ( + ;p^sm (— j , -2 = ^1 + . (5.6) 
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The next step is to write down the one-dimensional effective action for the closed string 
fields, namely the dilaton and the functions that enter the ansatz for the metric. This is 
similar to (2.9) but one has to take into account the ci-dependcncc of and that the DBI 
term is Nf times (5.3). As before, the WZ term does not contribute because it does not 
depend on the metric and the dilaton. We get: 



2b \ h b F J 

- \qIy2 - \Q)p{'yfe'^hS' -4g,cos^ (^) HiS.en^«F^) . (5.7) 

One can readily compute the equations of motion from this Lagrangian. An important point 
to note is that Xwv is not a field in this lagrangian. Its equation of motion (5.5) was derived 
from the worldvolume action. We find the following set of Euler-Lagrange equations: 

a^(iog6) = 0, 



h? 

dl{\ogS) = -26F^54 + 66F256-^g^e*6F2 56^os3^f=i + cos2^^ 

2 2 \:i2 2 

a^(logF) = AbF b -^e 65 p((T) — 2— e 06 cos — — sm — — , 

2 ^\ 2 2 

dl^ = Q)e''^bS^p{<Tf ^^Qfe^bS^'F'' cos' ^^{E2. (5.8) 
The first order constraint, generalizing (2.11), reads: 

~ ~2~^^ — ys — 2^ ^ 

^ m(i^^,i^^,iAn)_2,bF^s^+,bF^s^+ 

2b \ h S F J 

+ ^Q^4 + iQ/e'*^'^'M<^)' + 4Q/e*^''5'^'cos^^^ . (5.9) 
2 tl 2 2 ^2 

The only ingredient left to know is the precise expression for p{a). The answer is given in 
appendix C: 

p(^)=cos^^. (5.10) 

The system of equations (5.5), (5.8), (5.9), together with (5.10) defines the solution for a 
finite temperature D3-D7 background where the quarks are massive (all of them with the 
same modulus of the mass). If at some value a = Oq above the horizon Xwv reaches the 
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value Xwv — TT, it means that the flavor branes only have support for a > Gq. In that 
case, it has to be understood that equations (5.5), (5.8), (5.9), (5.10) are valid for a > aq 
whereas for o" < dg, one has the unflavored system of equations (namely one should substitute 
= for (7 < (jg in all the equations). At o" = (T^, one has to impose appropriate matching 
conditions. Detailed discussions of this feature in various supersymmetric cases can be found 
in [19, 20, 22, 23]. 

By deriving the constraint (5.9) with respect to a and inserting (5.5), (5.8), (5.10), one 
obtains an identity, proving that the system of equations is self-consistent. This can be 
considered a crosscheck for the whole procedure. The massless limit is given by Xwv — 0, 
Si = S2 = 1, and one recovers the set-up of section 2. 



5.1 The supersymmetric case 

We now consider the zero temperature solution of the system with massive flavors. One can 
check that the following first order equations solve (5.5), (5.8), (5.9): 

6 = 1, d,h^-Q,, d,F ^S'F (3-2^-^6" cos' ^ 

d,S = S'F' , d^xwv = -2S' tan ^ , 9.$ = S"e^' cos' ^ . (5.11) 

These are BPS conditions, generalizations of (2.12), which determine the supersymmetric 
solutions. It is convenient to define a p coordinate as dp = S'dcr. In the p coordinate, the 
embedding is just sin = e'''"'', where pq is an integration constant related to the bare 
quark masses. We can write explicit expressions for S, F, $, for p > pq?^ 

1 

S = a'^ e" ("l + e,(- + - p - -e^P-^'^P - -b'^P'^-^p + -e'p^-'P - -q'p^-'p* + e^^'-^^*)^ ' 
V 6 6 2 4 4 J 



,1 ^ il + e,{p,- p-e''P'^-''P + \e'P'^-'P + e^p''-''p* -le'P'^-'P*))'^ 

F — a'i eP Y 

(1 + e*(| + p^- p - ^e^Pi-^P - ^e^Pi-'ip + ^e^p<,-4^p - le^Pq-'ip* + g^^^-^P*)) ^ 

$ = $^ - log(l + e,{p,-p- e^P^-^P + ^e'P^-'P + e^P"-''"* - ^e'""-'"')) . (5.12) 

Setting Pq — > —00 one recovers the massless solution (2.13). 

We still have to write the solution for p < pq. The flavor branes do not reach this region 
and the equations of motion for the background are (5.11) with Qf = 0. Thus, the dilaton is 
constant and, by continuity, it has the value that can be read from (5.12) inserting p — pq. 

^jn = = log(l + (p, - p, - ^ + e^P^-^'^* - ^p^-^p*)) . (5.13) 



^^In integrating S,F,^, we have fixed integration constants as in the massless case. Moreover we have 
required S\p=p^ = -F|p=p, to ensure regularity in the IR (at p = —00) as explained in [20]. 
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The functions -S" and F are equal: 

S = F = a'"^ e^e-^(*^«-**) , {p < Pg) . (5.14) 

We defer a more detailed analysis of these solutions and of their generahzations to finite 
temperature to future work. 

6 Summary 

In this paper we have presented a backreacted supergravity solution dual to a flavored version 
of the Af = 4 SYM plasma, with massless flavors in the group U{l)^f. Our construction 
is also valid for black holes on AdS^ x X5 backgrounds, with X5 a generic Sasaki-Einstein 
space. In all the cases the T = solutions preserve A/" = 1 supersymmetry in 4d. 

The solutions we provide are perturbative, like the one presented in [25] for the unflavored 
thermal conifold theory. The relevant expansion parameter here, e/j ~ XhNf/Nc, weighs the 
internal flavor loop contributions to the unflavored field theory background. In the range 
of physical parameters associated to RHIC, we can reasonably expect values of €h ~ 0.24. 
On the other hand, this makes it a nice expansion parameter, and keeping corrections up 
to order ef^ yields a very good approximation to the full result. This is precisely the order 
at which we are able to provide a solution in the present paper. Moreover, leading order eh 
corrections are already of about 20% which are clearly dominant over curvature corrections. 
These come with a A^^''^ for string corrections to the supergravity lagrangian and with 
eh^h^ for corrections coming from the DBI action of the D7 branes [42]. For a realistic 
value Xh ~ 67r both of them are ~ 1%, although the precise value will very much depend 
on the prefactor, which can vary from one quantity to another. For example, for the free 
energy, the prefactor is of 0{1) [43], and therefore we expect the previous assertion about 
the dominance of flavor corrections to hold true for equilibrium quantities, like the entropy. 
On the other hand, dynamical quantities, like transport coefficients don't necessarily exhibit 
this hierarchy. Such is the case, for example, of the quotient ri/s which is uncorrected by the 
flavor, whereas curvature corrections give a large prefactor to X~^^^ of 0(10) [44]. 

Flavor corrected thermodynamics shows a departure from conformality at C(e^), which 
agrees with expectations, since we are adding massless flavors and therefore conformal sym- 
metry is broken by quantum effects. Being very mild, it would be very interesting to study 
the hydrodynamics of the plasma, in particular the bulk viscosity. 

We also considered the energy loss of probes in the plasma, by calculating the jet quenching 
parameter and the friction parameter. The outcome is that the energy loss is increased by 
the presence of matter in the fundamental representation. Even if this statement needs to 
be supplemented with some prescription of how to compare different theories, we found it 
to be true in all the reasonable comparison schemes that we analyzed. Thus, for example. 
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the value of the jet quenching parameter, as compared to the one in J\f — A SYM, is shifted 
towards the RHIC phenomenological window. It would be important to understand better 
the physical reasons behind this result and its implications, if any, for phenomenology. For 
example, there are indications that the smaller perturbative cross section of quarks with 
respect to gluons in the jet quenching process is not consistent with the data, suggesting 
that the strong coupling cross sections should be quite a bit larger (see for example [45], 
section 3.2). 

Besides the straightforward application to the study of phenomena like the ones mentioned 
above, having analytical control over a perturbative solution may prove instrumental when 
trying to obtain an exact backreacted one by numerical methods, so as to provide correct 
boundary and matching conditions. 

We finally began to attack the problem of the backrcaction of massive quarks. We derived 
the equations of motion and briefly considered the zero temperature case. It is conceivable 
that a solution at finite temperature can be calculated and analyzed following the methods 
in this paper. 

Apart from the ones outlined above, there are a number of properties and extensions of 
our solutions that deserve to be studied. Basically any observable calculated in jV = 4 SYM 
could be analyzed in the present setting. Some studies in this direction are currently being 
performed [46]. Also, one might think about putting the system at finite charge density and 
studying its properties. 
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A Technical details regarding the DBI contribution to 
the effective action 

Let us justify the form of the last term in (2.9), which comes from the DBI action of the 
flavor branes. In the supersymmetric case {b — 1), we can use the results of [7], where a 
two-form Q, related to the distribution density of the flavor brancs. was introduced. This 
form Q is such that dF^i) = —Qs^ = —gs'^i^^^K with the fi'-'-' being decomposable as 
the wedge product of two one-forms. Using supersymmetry, it was pointed out in [7] that 

Sdbi = -T7 / rfiOxv/=^e* E . ^J\n^^^n^^Qg^^gPQ. Since Vl can be directly computed 
from the expression of in (2.7), the T = version (6 = 1) of the last term in (2.9) is 
readily found. 

In order to generalize this result to the finite temperature case, one essentially needs to 
prove that the massless embeddings of the supersymmetric case are still solutions of the 
finite temperature case. With this purpose in mind, let us consider a generic AdS^-ERxX^ 
background with metric: 



ds'^h-^/^ [-bde + dx^]+h}'^ 



— + r^dsj^E + r^{dT + AKEf 



(A.l) 



Let us write the Kahler-Einstein metric as: 



(A.2) 



1=1 



where the one forms (Xj depend on the other angular coordinates ^a, ^ = 1,2,3. The one- form 
Ake will depend on a and and will not have a component along da. 

Let us consider a D7-brane extended along t,Xi,r,^a- It is trivial to show that r = const 
is always a solution of the embedding equations, as in the T = case. Picking a — a{r), the 
effective DBI lagrangian is proportional to: 



C ^ r^H{a, Ca) + b{r)r^a'^ , 



(A.3) 



where H{a, ^a) is the square root of the determinant of the metric of the three-dimensional 
compact space defined as the a,T — const slice of X^. 

The massless embedding a — const = ckq is a solution of the second order equations 
following from (A.3) if: 

daH{a,Q\a=ao^O ■ (A.4) 
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This is the same condition one would find in the T = (i.e. 6=1) case. Following an 

analogous reasoning, one can check that a = ao, r = const is also a solution of the D7 
equations of motion for the backreacted geometries (both T = and T > 0). Using the 
isometrics of the internal manifold, one can generate the family of embeddings needed for 
the smearing (both T = and T > 0). Thus, the last term of (2.9) for T > is a trivial 
generahzation of the T = result quoted above. 

Moreover, the relation between Qf and Nf can be determined by looking at this last term 
of (2.9). The DBI on-shell action for Nf "massless" branes is just: 

Sdbi = -T^Nf J d'x e* = -T,NfVol{X-,)V,,s j dae^h , (A.5) 

where VollX^) is the volume of the constant r, constant a = Oc slice of X5 that the D7 
wraps, as discussed above. In (2.9), the constant in front of this expression was defined as 
_^Q _^Yj^!d^^YL± _ Equating and using 2k\qT'j = gg, we find the second relation in (2.8). It is 

worth noticing that the factor appearing in (2.8) is just the volume transverse to any 

flavor brane. 



B Computation of the ADM energy and the free en- 
ergy 



Let us derive the result (3.9). The ADM energy is given by: 



.22 



Eadm = — \-V\9tt\ [ d^x^/detg^{KT - Kq) 



(B.l) 



The eight-dimensional integral is taken over a constant time, constant radius hypersurface. 
The symbols Kt and Kq are the extrinsic curvatures of the eight-dimensional subspace within 
the nine-dimensional (constant time) space, at finite and zero temperature, respectively. 



dfj_{y/dcFg^n^) where is a normalized vector perpendicular 



/Qrr 



-Ji!. The 



They are defined as K 

to the surface. Since we consider constant r hypersurfaces, this is just r?/ 
metric on the eight-dimensional slice must be equal for the zero temperature and finite 
temperature solutions that are being compared. With the solution written in (2.22), this 
is the case when r — without any further rescaling of coordinates. Using (2.5), (2.15), 
(2.20), the explicit expression for the ADM energy reads: 



E 



ttN^V. 



Wol{X5)R» R 



1 - 



dr{r^R^S^F) 



RS^F 



^ ^ Or 



„3 p2 c4 17 ^ 



RS^Fo 



(B.2) 
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We use the notations in [47]. 
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Plugging in the solutions (2.17), (2.22) we get: 



1 + ^ + 0{el) 



(B.3) 



where we have discarded terms suppressed by powers of We can now write the result in 
terms of eh and T by substituting (3.2) and (3.4), and obtain (3.9). 

Let us now directly compute the free energy in (3.11). It is given by the Euclidean action 
I (renormalized by subtracting the zero temperature result) evaluated on the solution and 
divided by the inverse temperature /3 [48]. We can write: 



(B.4) 



We take as the radial cut-off for the integrals, such that the finite T and zero T geometries 
coincide. The only subtlety comes from the fact that gu{r^) ^ g^^{r^). Thus, in order to 
compare the solutions, we ought to rescale the Euclidean time of the zero temperature 

solution such that its period is [3q = /3 ^f** ^* = P-Jl — where /3 is the period of the 
Euclidean time of the finite T solution. 

The bulk action comes from the sum of (2.3) and (2.4). The WZ term in (2.4) does 
not contribute on-shell for the massless embeddings we are discussing here. By using the 
equations of motion of the system, one can prove that the integrand of the on-shell action is 
just equal to — 8r^. We get a very simple expression: 



"2/^10 



V^Vol{X^) 



— 



bulk,T — ->-bulk,T 



-o) 



r 



4 /■?■* 

h 



(3 1 {-8r')dr - (3Jl- ^ {-8r-')dr 



(B.5) 

In this expression and the following, it is understood that corrections in powers of ^ and of 
order or higher are discarded. The surface action comes from the Gibbons-Hawking term: 



■'surf 



K 



= ^ / KdT, , 



(B.6) 



10 -'S 



where E is the 9d subspace at r = r* and K the extrinsic curvature of this hypersurface 
within the full lOd space. Exphcitly: 



V^Vol{X^) 



{Xsurf,T — Xsurf,T=o) — 2 



r\/l - ^ 
RS^F 



(3 dr{r^R\ll - "^S^F) 



(B.7) 
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Summing the bulk and surface contributions and inserting (2.5), (2.15), we get: 



1 



In terms of and T from (3.2), (3.4), we recover (3.11). 



(B.8) 



C The smeared Karch-Katz model 

In this appendix, we want to write down the family of embeddings that participate in the 
smearing of massive flavors in AdS^ x (at the end of the appendix we comment on the 
generalization to finite temperature backreacted backgrounds). They are all holomorphic 
embeddings in a given set of complex coordinates and therefore the overall preserved super- 
symmetry is jV = 1. From the geometric point of view, rewriting the Af = 4 theory in A/" = 1 
components amounts to rewriting the metric of in the canonical Sasaki-Einstein form, i.e. 
as a U{1) bundle over a complex 4d Kahler manifold. In this case, the latter is CP^. 

In the following, we use the notation from [49]. Let us consider (D^ with complex coordi- 
nates Z^, Z'^, Z^ and metric: 

ds^ = + \dZ^\^ + \dZ'\^ . (C.l) 

We introduce a radial variable r such that = rz', with = and the spanning 

a unit five-sphere: 

+ i^T +kT = 1 ■ (C.2) 

This constraint is invariant under SU{3) rotations of the z^ and under the transformation 
_^ e'^'^z^, where a is a phase. The space CP^ is the space of orbits under the action of 
this U{1) group. We can parameterize: 

z = COS — COS - 62 '^v'-rv^; ^ ^ _ gjj^ _ ^2 ^ ' ^i" z = sm — e . 0.3 

2 2 ' 2 2 ' 2 ^ ^ 

By inserting (C.3) into (C.l), one gets the metric of as written in (2.1), (5.1). The gauge 
theory U{1)r symmetry is related to the isometry generated by the Killing vector 9^. The 
volume element of 5"^ in these coordinates is: 

£[S^^) = -L sin 6' sin x(l + cos x) dr A dip A d(p A d9 A dx ■ (C.4) 

The range of the angles is < X; ^ ^ tJ"; < r < 27r, < ip < in such that / £{S^) = 
VnI(S^^) = IT'''. AVe now discuss the holomorphic embeddings. We can take Z^ — Vq with Vg 



^^Oi course, rewriting Z'' = + iy^~^^ and so on, with real, the defining equation for reads 
V'' — 1- This is actually invariant under the larger group S'0(6) ~ SU{4). The choice of an holomor- 
phic parameterization breaks this symmetry group to the smaller SU{i) x U{1)r one. 
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real (and related to the modulus of the quark mass) as the "fiducial" (Karch-Katz) embedding 

[50]. In terms of the coordinates, the fiducial worldvolume satisfies sin = y and r = 0. 
Notice that r is constrained by r > and that it diverges for x = 0- Acting on it with the 
SU (3) X U{1)r isometry group, we can write the following family of holomorphic embeddings: 

3 3 

J^aiZ" = rqe'^ , with ^ |ai|2 = 1 , (C.5) 

i=l 1=1 

where < /3 < 27r is associated to the phase of the quark mass and the are complex 
parameters. They span a unit five-sphere and therefore we can represent them by using a 
parameterization similar to the one for the z's in (C.3), in terms of the tilded coordinates 
X, 9, f, ■0, (p. Then the generahzed Karch-Katz embedding (C.5) can be rewritten as: 



gi(r+f)gi(V'+V!+</'+'?)p — Ll^iP (C.6) 

r 



with: 



1 = COS — COS - COS — COS — h COS — sm - cos — sm - e -\- sm — sm — e 2 vf^-rv'-Tv^-rv^; _ 

2222 2222 22 

(C.7) 

The supcrpotential of the theory reads schematically: 

W = $i[*2, *3] + g(ai$i + a2$2 + oa^a - rn)q , (C.8) 

where $1,2,3 a-re the three complex scalars of jV = 4 and the g's are the matter multiplets of 
mass m ~ VqC^^ in the fundamental representation provided by the D7-branes. 

Following the procedure of [20], our goal now is to consider a homogeneous, symmetry 
preserving, distribution of D7-branes using this six-parameter family of embeddings and 
compute their charge density f2, such that dFi^i) = —Qs^- Let us start by splitting the 
complex equation (C.6) into two real ones /i = /2 = with: 

/i = 2(T + f)+z^ + z^ + (^ + (^ + 2Arg[r]-2/3 + 47rn, f^^\r\^--^. (C.9) 
The fl is built by summing the 2-forms locally orthonormal to each brane worldvolume [20]: 



n 



2Nf 

^ sin x(l -|- cos x) sin 6 



(47r) 



S{fi)S{f2)dfi A dxdedfdi^d^dp , (C.IO) 



where the factor in the square brackets corresponds to the symmetry preserving normalized 
density of the fiavor branes. In principle, performing the integral (C.IO) would be a challenge 
but, because of the preserved symmetries, we know that = Qfp{r){dT + Acp^), as in 
(5.2), and thus: 

1 Nf 

n ^ = --^ [2p(r) Jcp2 + (drp(r)) dr A (dr + Acp^))] . (C.ll) 
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Thus, the problem of solving the integral (C.IO) is reduced to finding a single function p{r). 
The easiest way to proceed is to evaluate from (C.IO) the dr A dr component, at x = tJ"- 
A straightforward analysis yields (see [20, 22, 23] for details of analogous computations in 
different cases): 

p(r) = , (r <rg) , 

p{r) = (l-^) , ir>rg) • (C.12) 

As expected, p{r) vanishes at scales below the quark masses while it asymptotes to 1 in the 
UV. 

Up to here, we have written the family of embeddings (C.6) and the expression for the 
corresponding charge density (C.ll), (C.12) when the background is not backreacted and 

the temperature vanishes. However, the computation is easily generalized to cases with 
backrcaction, finite temperature, or both. This is because the isometrics and consequently 
the integrals over the angles arc exactly the same. One should just change slightly the 
expression for /2 in (C.9) to be /2 = |r|^ — sin^ x^^(cr), where Xwv{<^) parameterizes the 
fiducial embedding and, as in the main text, we have used a for the radial coordinate. The 
argument that led to (C.12) is readily generahzed to obtain (5.10). 
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